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Abstract 

We predict a gigantically long lifetime of the first excited state of an interstitial lithium donor 
in silicon. The nature of this effect roots in the anomalous level structure of the Is Li manifold 
under external stress. Namely, the coupling between the lowest two states of the opposite parity 
is very weak and occurs via intervalley phonon transitions only. We propose to use these states 
under the controlled ac and dc stress to process quantum information. We find an unusual form 
of the elastic-dipole interaction between different donors. This interaction scales with the inter- 
donor distance R as R~'^ or R~^ for the transitions between the states of the same or opposite 
parity, respectively. The long-range R~^ interaction provides a high fidelity mechanism for 2-qubit 
operations. 
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A broad effort is currently underway to develop new techniques and systems for quantum 
information processing (QIP) that involves coherent manipulation and control of a large 
array of 2-level systems (qubits) for quantum storage and computing purposes. Implemen- 
tation of the quantum computer (QC) as a solid-state device has a number of unique advan- 
tages related to scalability and the prospect of integration with state-of-the-art semiconduc- 
tor device fabrication technology. One of the main challenges in developing a solid-state QC 
is the requirement of a relatively low qubit decoherence rate. The natural candidates in this 
respect are QC schemes based on spins as quantum bits, e.g. nuclear spins of '^^P donors ll, 
or electron spins in semiconductor quantum dots [3]. While spins are relatively well isolated 
from the environment and characterized by long decoherence times, the implementation of 
the spin-based QC is complicated by a short-range character of exchange interaction. This 
provides additional difficulties to very challenging problems of a single spin measurement and 
control as well as nm-precision fabrication of a periodic qubit array. Additionally, the clock 
frequency for nuclear spin QC is relatively low. The optically enhanced RKKY mechanism 
for electron spin-exchange p| in semiconductor quantum dots allows for greater inter-qubit 
distances, but it is not clear at present how to make this QC scheme scalable. 

An interesting alternative to the spin-based qubits is provided by charge-based qubits 
formed by the two lowest states of a single-charged pair of quantum dots jj, ^ , or donor 
atoms 0] . Those systems do not depend on the single-spin readout and rely on a long-range 
electric-dipole interaction between the qubits for performing quantum logic operations. The 
QC scheme proposed in uses acceptors in silicon to encode qubits with a decoherence 
rate ~ 1 KHz and relies on a long-range elastic-dipole interaction between qubits (~ R~^). 
A serious source of decoherence in all charge-based qubits is due to charge fluctuations in 
the surrounding environment ^. In particular, relaxation of the charge traps can produce 
a shot noise signal leading to decoherence times that are significantly shorter than that in 
the spin-based QC sl- 

In this Letter we exploit the multivalley nature of the Si conduction band and anomalous 
structure of the ground state of an interstitial Li donor to demonstrate that under certain 
conditions the first excited state of the single Li donor can have a lifetime comparable to 
that of a nuclear spin (in access of 1 sec). The two lowest Is states of the Li donor in Si 
can be used to encode the qubit that is extremely well isolated from the electromagnetic 
fiuctuations of the surrounding. 
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The conduction band of Si has six valley minima kj located along six equivalent directions, 
[x, —X, y, —y, z, —z] at about 85% of the distance 27r/asi to the Brillouin zone boundary (here 
asi is a lattice constant for Si). In the framework of the effective mass theory (EMT) the 
ground state of a shallow donor electron is 6-fold degenerate. The valley-orbit interaction 
removes the degeneracy and gives rise to the splitting of the donor levels Q . Excited states 
within the ground-state manifold are produced by splitting of the Is hydrogenic energy 
levels and decay predominantly via emission of acoustic phonons. The donor electron wave 
function can be written in the standard form = X]fj|=i '^j Xlk g ^I+gI'^)- Here k is a 
wave vector in the first Brillouin zone, G is a reciprocal lattice vector, |k) are the Bloch 
functions and index fi labels the irreducible representations of the point group that are 
characterized by the coefficients aj. Also, are the hydrogenic envelope functions in 
/c— space that are strongly localized in the vicinity of the corresponding valley center kj and 
the summation above is taken over all valleys j = ±1, ±2, ±3. 

We now consider the acoustic phonon driven transition rate from a state l^^'*') to 
Wf,f,: = ^TT/h'^Y.qu ~ ^q'^)' where hw^^^f, is the change in the donor energy 

during the transition, Qq^, is the frequency of an acoustic phonon mode u with the wavevector 
q, and V^^'(qi^) are the matrix elements of the electron-phonon interaction Hamiltonian, 

(^'^li/el-phl*'^') = Eq. l^MM'(q^)(V + where 
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K''= E «r<T^k*^u+G-Mkq.. (1) 

|i|,|j|=l k,G 

Here Aikqu is the matrix element of i?ei-ph between the Bloch functions |k) and |k + q) [lol |. 

For the long- wavelength acoustic phonons, with g <^ |kj — kj|, the dominant contribution 
to (0) comes from the intravalley terms with i = j. The only other important terms 
correspond to the intervalley umklapp processes between the valleys j and —j centered at the 
band minima kj and k'_^- = h j + Gj, where k_j = — k^, = {An / asi)kj, and kj = 'kj/kj. 
The points kj and k'_j lie on the opposite sides of the Brillouin zone boundary and correspond 
to the closest possible intervalley separation kq = O.Gir/asi- All the other intervalley {i ^ j) 
and umklapp terms (G 7^ 0) are much smaller, since they are proportional to the overlap 
of the envelope functions from much farther separated valleys. For both the intravalley and 
intervalley {j —j) umklapp processes the sum over k in Eq. is dominated by the small 
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vicinity of the valley centers, where A^kq;^ ~ ■M±kjciu = M^^^, and 

Ml, = {n/2pVn^,f^ ■ q)(% ■ e^,) + Hd(q ■ Cq,)] . (2) 

Here "E^ and are the deformation potential constants, ec^u is the polarization vector of 
a phonon mode qz/, p and are the mass density and volume of the crystal, respectively. 
Finally, we obtain 
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KT' = EE "f^k [«f + «^^>kiq+Gj , (3) 

b1=i k 

The EMT donor eigenf unctions \l/'^(r) have a certain parity which is equal to the parity 
of the multiples Q^j^k under the operation kj For example, the parity is even for 

the s-type singlet or doublet donor states belonging to the irreducible representations Ai 
or but it is odd for s-type triplet states belonging to T2 etc. In the limit of small q the 
intervalley processes described by the second term in Eq. (jS)) give rather small contribution 
and presumably can be neglected. Then, as it follows from Eq. ((2)) and Q, the matrix 
elements between the states of the opposite parity, V^l^' = 0, due to the cancellation of the 
intravalley j and —j terms. Based on this observation one can expect that the probabilities of 
the even-odd transitions involving long-wave acoustic phonons are much smaller than those 
of the same parity transitions. Analogous symmetry arguments explain a strong suppression 
of the ls{Ai) — > \s{T2) Raman transition in silicon donors One can also show that 

if the intervalley processes in are neglected the transitions between the opposite-parity 
donor states are forbidden in all orders in V^^^' . 

We now consider implications of this effect for the lithium donor in silicon. Li is an 
interstitial impurity with the T^-site symmetry. It has an anomalous fivefold degenerate 

lies above, with the splitting 



Ij] . If a uniaxial compressive 



ls{E + T2) ground state while the fully symmetric state 
between the states conventionally denoted as 6Ac =1.76 meV 
stress Fz is applied along the (001) direction the site symmetry is reduced from to 
and has a distinct mirror-rotation axis /cs. Then the ground state level is split into three 
levels as shown in FigH jl^ . Q]- The new ground state denoted as |0) has an odd parity 
with the coefficients = (0,0,0,0, 1, — 1)/-\/2. The first excited state, |1), is a singlet, it 
has an even parity, aj^"* = (b, b , b, b, a, a) where a= a{e), b=b{e). At small magnitude of the 
applied stress \Fz\ one has 13|: a = (6 + e)/6-\/3, 6 = (3 — e)/6\^, where 

e = Euisn - Si2)F,/3A„ (4) 



and Sii, Si2 are elastic compliance constants. The second excited energy level is a triplet, 
it consists of an even parity state, |2), with aj^'* = 1, —1, —1,0,0) and two odd parity 
states. 

From the above it follows that only the second term in Eq. ^ with i = —j = ±3 
contributes to the phonon decay rate Wio of the state The decay rate calculated with 
the Kohn-Luttinger envelope functions joj reads 



Wio = r,o - u , i^.... /o^216 • ^u...7 ■ (5) 



35 [1 + (a||Ko/2)2]6 TThfm 



Here Ut is the transverse speed of sound, ay is the longitudinal Bohr radius for the is donor 
state in Si, and we assumed that 2T\:ut/uJu) ^ a\\. We emphasize a very fast decrease of 
the transition rate between the opposite parity states with the level separation, Wiq oc tUgi. 
This is in contrast with the conventional dependence of the transition rate between the same 



parity states 
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IJ]. For instance, W21 oc UJ21 for the transition |2) We note that the 

decay rate oc uj^ was found in jl^ for the electron transition in a singly- charged double-dot 
system between the antisymmetric excited state and symmetric ground state separated by 
the energy gap Ituj. However, a key distinction of the expression Q from the results of 
Ref. is that it contains an additional small factor, 8a^(£:)/35(a||fi:o/2)^'^ ~ 10"^, due to 
the overlap of the envelope functions from different valleys. It is the combination of these 
two small factors in (jSJ that gives rise to the giant lifetime of the state |1) at the low stress 
values (see Fig. [U). This is one of the central results of this work. 



In addition to the energy decay there also exist the dephasing processes |7( due to quasi- 
elastic phonon scattering off donors. Detailed analysis shows that these are not important 
at low temperatures. For example, for i^ei-ph given by Eq. the dephasing rate Wio = 
uolT/ToY^, with i/q ~ 2 ■ 10^^ Hz and Tq = hut/a±=19 K. Thus at T=0.1 K the dephasing 
rate is negligibly small. Finally, the dipole moment of |0) transition is less then 10"'^ eA 
and therefore its quantum phase is very robust against the charge noise. 

The fact that electronic states of a Si:Li donor are extremely well-isolated from the envi- 
ronment provides an interesting possibility to consider this system for quantum information 
processing with a pair of states |0), |1) used to encode a qubit. In what follows we make 
use of the fact that by varying the local pressure Fz on a given qubit one can change the 
frequencies u^^/ and therefore lifetimes r^^/ of the corresponding transitions. As seen in 
Fig. ^ there is a dramatic difference between the lifetimes of the even-odd and same-parity 
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FIG. 1: Dependence of the lifetimes rio, T21 of the first two excited states of Li donor in Si on 
e. The insert shows the energy levels of an isolated Li donor in Si under the compressive force 
F||(001). Bohr frequencies wio = a;2i/2 = eAc//j. 

transitions which exceeds seven orders of magnitude at e = 0.5. For smaller e the lifetime 
Tio becomes extremely large, it reaches seconds and even minutes for e < 0.2 or fuvm < 0.06 
meV. This energy scale can be used for a qubit level separation. 

One-qubit operations that are required for quantum computation with Si:Li donors can 
be implemented by applying tailored pulses of stress or electric field locally to selected 
donor atoms (Fig. E)). For example, a rectangular stress pulse of duration t^c in (001) 
direction will shift donor transition frequency loiq and produce a phase gate exp[iT^c 10 Sz]- 
Here ^2 is a ^-component of the pseudo-spin- 1 operator S acting upon the qubit states 
|0) = II) and |1) = If). One can also apply a low-amplitude periodically modulated stress 
pulse v4^ooi> coscJiot of duration Tac in (001) direction. This pulse will resonantly drive the 
transition |0) — |1) of a given qubit. It will produce a pseudo-spin rotation around x-axis, 
X{(p) = exp[i(pSx/2], ip = 2tiQx- An estimate for fl^ in the rotating wave approximation 
gives 

ilx = -= + r.d)KQa\\. (6) 



6 



Here ui is the longitudinal speed of sound). With a modest ac stress amplitude A(ooi) — 
10^ dyn/cm^ and a driving frequency u;io/27r=10 GHz one gets fix-/27r ~630 MHz. Then 
the TT-pulse X^ir) will have a duration ti ~ 0.4 ns and a corresponding decoherence-induced 
error ~ tiW^iq will be less then 10~^. 

Another controlling option consists in applying pulses of the time- dependent electric field 
E(t)||(001) caused by the modulation of a bias voltage between the local electrodes {A,B 
in Fig. 12)). The electric field will produce a quadratic Stark shift of uio and will induce the 
Rabi oscillations between the qubit states |0) and |1) via the virtual transitions through the 
2pQ state. For an electric field amplitude E ~ 10'^ V/cm the period of the Rabi oscillations 
Ti oc E""^ is of the order of 0.1 ns. 

Li donors in Si are coupled to each other via acoustic phonon field and behave as elastic 
dipoles. As in the case of the shallow acceptors in Si 7|, this coupling can be used as a 
"data bus" to implement two-qubit logic gates. The matrix elements of the phonon-mediated 
coupling Hamiltonian, Hij = ^ ' a a' (K/^fl + between a pair of donors at the 

sites Hi and Rj can be expressed as follows: 

where we assumed that ct;^^/|Ri — Rjl/wj ^ 1. 

Consider first a resonant excitation transfer (RET) between donors i and j involving a 
pair of transitions, /i — /i' on site i and fi' —* fi on site j, that are in resonance with each 
other. The coupling constant corresponding to this process hgl^l^ = G^^, For the RET 
involving the even-odd transitions |0) ^ |1) this coupling constant reads 



^ = 3 - 4 ■ 4a + 5 ■ — . 8) 



Here a = S^^/S^ and we assumed that Rjj = Rj — Rj lies in the plane normal to (001). We 
note an unusually steep falloff of the coupling constant with the donor-donor separation. 
In fact, one can show that if the RET involves transitions between the states of the same 

coupling constant for the |1) ^ |2) transition in Si:Li is 

g■}/W2l = liut/u;2^R^,f (9) 

where 7 = (5/16) (2 + 7(1 — uf/uf)) ~2.0. This long-range dipole interaction is due to the 
intravalley terms in V^^ that are absent in V"q°. It can be used to selectively generate an 



entanglement between two donors by tuning their energy levels into resonance with each 
other. 

Consider now a diagonal part of the phonon-mediated donor-donor Hamiltonian generated 
by the diagonal matrix elements V^l^ in Eq.((7)). This is an Ising-type dipole-dipole interaction 
resulting from the dependence of a polaron shift in the electron energy of a given donor on the 
electron states of the neighboring donors. The effective interaction Hamiltonian projected 
onto the states /i = 0, 1 of donors located at Rj has the form: 

^(^) = 2 5Z ['^^^■^- + 9^MS- + h.c] + H,{t). (10) 



Here S'^ and are raising and lowering pseudo-spin operators for jth qubit, Hj{t) is 



a single-qubit Hamiltonian, gij = gjj, and the Ising-exchange coupling constant hJij = 
11 + ^00 00 ~ 2(^00 1]^. Due to the symmetry of the problem, Jij does not depend on 
the in-plane orientation of Rjj_L(001). Using (0), ©j© and the form of the valley-orbit 
coefficients one can obtain 



^'^-32nhfm^^ [ ^ + 3 ' «f J ' ^^^^ 

Here the 1 / Rfj law is also due to the intravalley terms that are always present in the diagonal 
matrix elements V^Jf". For Rij = lOOnm, Jij/ir ~ 10 MHz while gij/n=OA Hz. Because of the 
detuning of the donor transitions from resonance the RET frequency at |0) |1) transition 
is further reduced and can always be neglected. Contrary to that, the dipole interaction in 
Ti.(t) always persists. It is a basic long-range 2-qubit interaction in the liquid-state NMR 
and also in the solid-state QC scheme based on magnetic-dipole coupling of donor spins in 
silicon llSll . The Hamiltonian (fTUI) with g^j = is sufficient to execute any one- or two-qubit 
gates 15[- It serves as a basis of our QC scheme. 

The QC scheme utilizing Ising-type dipole coupling has been described in Ref. It is 
based on the refocusing technique where an auxiliary set of vr pulses is applied to a group 
of qubits j = 1 : n. The vr-pulses cancel the dipole interaction terms in pO|l for all members 
of this group except for a selected pair of neighboring qubits. The error of two-qubit gates, 
caused by the non-exact cancellation, p ~ {ti/t2)^ |15[, where ti and T2 = ir/Jjj^i are 
characteristic time scales for one- and two-qubit gates, respectively. For Rjj+i = 100 nm 
T2 = 10~^s and ti <10 ns one has p < 10^^. The error of two-qubit gates due to the 
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phonon-decay of the state |1) is extremely small. The corresponding quality factor q for the 
two-qubit gate 



is probably greater then 10^ at a;io=0.06 meV//i. 

We note that a frequency of resonant ac stress modulation oj ^ uiq is typically limited 
by a high-frequency cutoff of the response function of a nano-scale piezoelectric transducer 
(cf. Fig. 12)). One can show that for practically feasible cases this implies uoiq < 10^'^ rad/s. 
However one could keep the qubit energy difference at much higher levels, ujiq ~ lO^^rad/s 
(0.1-0.06 meV), at all times, except the short time intervals when a given qubit is involved in 
a gate operation. Then the energy difference can be adiabatically reduced the smaller value 
ujiQ < 10^° rad/s. If the time between gates involving a given qubit is ^ W{f^ the qubit 
energy levels will not be thermally repopulated and the working temperature of a quantum 
computer can be set to T* ~ fWJio/kB — 50-100 mK. 

Similarly, one can use the adiabatic stress reduction in the alternative three-level QC 
scheme where a long-range RET, based on the donor transition 1 — 2 (see above) with 
Rabi frequency gfj (Eq. Q), is used to execute 2-qubit gates. This process provides the 
most direct tool to generate entangled states c|ljOj+i) + c'|Ojlj+i) between a given pair 
of neighboring donors {j,j + 1) whose energy levels are brought into to resonance using 
locally applied stress pulses. To achieve high-fidelity 2-qubit gates within this scheme we 
need much smaller interlevel separation hu!2i than hu!2i =0.12 meV at e =0.2. Therefore 
whenever we need to execute a 2-qubit gate we will adiabatically reduce the stress e{t) from 
Eq = 0.2 to 82 — 0.002 {Fz ~ 1.3 ■ 10^ dyn/cm^) which will set the transition frequency to 
cj2i/27r ~ 0.24 GHz. 

Within the three-level scheme the Rabi frequency gf]_^.i defines the QC clock frequency. 
In analogy to the above we define the quality factor q as the ratio of gf]^i to the phonon 
decay rate of the state |2) at the finite temperature T*: 



where R = -Rj,i+i and n2i(T*) = [exp(}TiU2i/kBT*) — 1]~ is the Planck's distribution func- 
tion. Consider a particular example with R=50 nm, e = 2 ■ 10""^, and T* =100 mK. In this 
case the characteristic time for a 2-qubit gate, 27r/5'i2 ~0.2 fisec, is much shorter than the 



q = 



T2W10 



1 



35 [1 + (a||/to/2)Y Ut (-1 + 5/3 jut/uif) 



(12) 




(13) 
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FIG. 2: Quality factor for the 3-level quantum computing scheme: dependence of the operational 
temperature T* on the inter-qubit distance R for different values of the quality factor q and 
= 1.3 • lO^dyn/cm^ {hwu = 0.001 meV) 

relaxation time T21 ^ 3 msec and therefore the donor electron levels will not be thermally 
repopulated during RET despite ksT* /fujji2 ~ 8. As soon as the two-qubit gate is executed 
the qubits will be adiabatically returned into their high-stress state. With this choice of 
parameters the QC clock frequency is (72i/27r ~ 5.2 MHz and the quality factor of a 2-qubit 
gate 5'2i/W^2i ~ 10^. Solving Eq. (fr5|) for T* we obtain the operational temperature as a 
function of desired quahty factor q and inter-qubit separation R: 

T* = -W3fcBlnfl-f -^^y (14) 

V 37 y 

By lowering the operational temperature T* < 100 niK we can significantly increase the 
quality factor or/and increase the interqubit separation R as shown in Fig. |21 

By varying the locally applied stress and voltage on the local electrodes (see Fig. Ej) 
one can selectively shift the \s{E + T2) — 2po transition frequencies of a given donor atom 
by h^uj ~ 0.3 meV without significantly affecting the qubit level separation fvjJiQ. Then the 
donor excitation lines ls(-E' + T2) — 2po can be selectively brought into resonance with a 
globally applied infrared radiation field since their homogeneous broadening (~55 MHz [ifij l) 
is ~500 times smaller then the frequency shift Au;. Such a selective resonant excitation can 
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be used to perform a measurement of the qubit state. Indeed, according to the polarization 
selection rules 3] a linearly polarized electric field E±F will resonantly excite the transitions 
from the state |1) to the states of the 2pQ manifold but will not affect the donor electron 
that was in the state |0) before the IR pulse arrived. Detection of the dipole moment of this 
state (~ 40 Ae) by means of a single electron transistor (SET) measurement will determine 
whether the state |1) was occupied prior to the excitation. The fabrication of the Si:Li QC 
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FIG. 3: Schematic of the Si:Li-based quantum computer. An etched piezoceramic structure (blue), 
which is placed on the top Si layer (green), has a row of columns of a width w ~ 20 nm separated 
by arches. The columns are used to generate a local (001) pressure, controlled by the voltage biases 
between the pairs of B-C electrodes. The decay factor of the local stress (and electric) field over 
the interqubit distance d ~ 100 nm is ~ {w/2df = 10-2 . The residual stress and electric field 
at each Li donor arising from the local control of the neighboring qubits can be compensated by 
the synchronized adjustment of voltages for all qubits. 



will utilize very large electro-mobilities of Li ions in Si (at T=400K and i?=10^V/cm the 
drift velocity is 0.4 cm/s). By controlling the electric field between the electrodes A and B 
(Fig. ini) the Li ions, sequentially deposited on the surface of a super-pure Si layer (positions 
1), will be dragged to their target positions 2 to form a subsurface array with the period 
d ~ lOOnm. Subsequent lowering of the temperature to its operational value T ~ 50-100 
mK will form a stable regular array of the neutral Li donors. 

In conclusion, we have shown that the valley-orbit structure of the ground and first excited 
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states of the Si: Li donor gives rise to a gigantic lifetime of the excited state. We have found 
that the elastic dipole coupling between two donors scales with the inter-donor distance R 
as either \/R^ or depending on their electronic states. We build on these findings 

and propose a QC scheme with the spinless stress-defined elastically coupled qubits. It is 
characterized by extremely small error rate (lesser then 10^^), relatively large operational 
temperatures and interqubit separations that are 5 to 10 times greater then those used in 
the nuclear spin solid state QC. 

We acknowledge valuable discussions with M. Dykman, B. Golding, and M. Foygel. This 
work was supported in part by NASA Revolutionary Computing Program (V. S.) and NSF 
(A. P.). 
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